Displacement-enhanced continuous-variable entanglement concentration 
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We study entanglement concentration of continuous variable Gaussian states by local photon 
subtractions enhanced by coherent displacements. Instead of the previously considered symmetric 
two-mode squeezed vacuum states, we investigate the protocol for input states in the form of split 
single-mode squeezed vacuum, i.e., states obtained by mixing a single-mode squeezed vacuum with 
a vacuum state on a beam splitter, which is an experimentally highly relevant configuration. We 
analyze two scenarios in which the displacement-enhanced photon subtraction is performed either 
only on one, or on both of the modes and show that local displacements can lead to improved 
performance of the concentration protocol. 

PACS numbers: 03.67.Bg, 03.67.Hk, 42.50.Dv 



I. INTRODUCTION 

In recent years, a great deal of effort has been dedi- 
cated to the study of continuous variable quantum infor- 
mation processing (CV QIP) that instead of qubits uti- 
lizes infinite dimensional Hilbert spaces of field modes. 
Of particular importance is the class of Gaussian states 
and Gaussian operations P, 0| which is experimentally 
feasible by linear optics, parametric amplifiers, and ho- 
modyne detection. Several important protocols such as 
quantum teleportation and quantum cryptographv can 
be implemented using only Gaussian means |3|45| and 
have already been demonstrated experimentally [a 0| ■ 

One of the key advantages of CV QIP is that the en- 
tangled quantum states ^ , which represent a central in- 
gredient of many quantum information processing tasks, 
can be generated deterministically using optical paramet- 
ric amplifiers and beam splitters. However, entanglement 
distillation and concentration of Gaussian states belongs 
to the tasks that require some non-Gaussian operation 
[9l-[Tl| . Recall that entanglement distillation and concen- 
tration is a crucial tool for entanglement enhancement 
and suppression of losses and noise in quantum commu- 
nication. This procedure enables two distant parties to 
extract, from a large number of shared weakly entangled 
states, a smaller number of highly entangled states by 
means of local operations and classical communication 
p"3 - [T^ |. Examples of non-Gaussian operations suitable 
for concentration of Gaussian entanglement include pho- 
ton subtraction [isl - fisj , photon addition [23| , inter- 
ference with ancillary single photons |2l| - |23j . and Kerr 
nonlinearity [13, [11] . 

Very recently, it was shown that, despite the fact 
that Gaussian operations cannot lead to concentration of 
Gaussian entanglement, they can improve performance of 
entanglement concentration based on photon subtraction 
[26r,28j . Here, we follow this line of research but instead 
of a symmetric two-mode squeezed vacuum state, consid- 
ered as an input state in Refs. [2^, [2^, we analyze the 
protocol for a split single-mode squeezed vacuum state, 
i.e., single-mode squeezed vacuum mixed on a beam split- 
ter with a vacuum state. This latter configuration is 



highly relevant because the recent experimental demon- 
stration of entanglement concentration by photon sub- 
traction was performed precisely with this kind of entan- 
gled state [231 • We assume a combination of local photon 
subtraction with coherent displacements, since squeezing, 
proposed in Ref . (26j , is experimentally much more chal- 
lenging. Following the two entanglement-concentration 
schemes implemented by Takahashi et al. [29t , we inves- 
tigate configurations in which the photon subtraction is 
performed either by only one of the parties, or on both 
modes simultaneously. We show that while in the former 
scenario the displacements do not improve entanglement 
concentration, in the latter case the output entanglement 
is enhanced by the use of local displacements. 

The paper is organized as follows: We introduce the 
entanglement concentration protocol and the considered 
input states in Section [ll| The main results are presented 
in Sections IIIII and ITVl where the performance of the pro- 
tocol is studied, first for weakly squeezed input states and 
later for states with arbitrary squeezing. Optimal dis- 
placements in dependence on the initial squeezing are de- 
termined, together with the resulting entanglement and 
success probability of the protocol. Finally, we conclude 
in Section fVl 



II. ENTANGLEMENT CONCENTRATION 
PROTOCOL 

The considered entanglement-concentration scheme is 
illustrated in Fig. [TJ The input entangled Gaussian state 
is generated by splitting a single-mode squeezed vacuum 
(SMSV) state into two modes on a beam splitter BS. To 
simplify the subsequent analysis and obtain maximum 
insight into the studied protocol, we consider a pure input 
squeezed vacuum state with Fock-state expansion. 



V cosh s 



2"n! 



(tanhs)"|2n), (1) 



where s denotes the squeezing constant. After mixing 
with vacuum on beam splitter BS with amplitude trans- 



mittance t and reflectance r, we obtain the two-mode 
state 
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where A = tanhs and \m,n) is a shorthand notation for 
\m)A ® \n)B- 

The two modes A and B are then distributed to two 
parties, Alice and Bob, who attempt to increase the en- 
tanglement of the shared state by local coherent displace- 
ments Zj^i (a) [Z)b(/3)] and single photon subtractions. As 
schematically indicated in Fig. [1] photon subtraction is 
achieved using a highly unbalanced beam splitter which 
reflects a small part of the incoming light onto a detector 
whose click heralds successful subtraction [30l - [3^ . Here 
and in Section IIIIl we model photon subtraction by the 
action of an annihilation operator d (b) onto the state. A 
more realistic treatment taking into account reflectance 
of the tap-off beam splitters and properties of the detec- 
tors will be provided in Section IIVI It will prove help- 
ful to assume that, after photon subtractions, Alice and 
Bob undo the displacements by applying inverse displace- 
ments D\{a) = Da{—oi) [Dg{(3)] to the shared state. 
Note that these latter local unitary operations do not 
change the amount of entanglement so they are not a 
necessary part of the entanglement concentration proto- 
col and can be skipped in practice. 

In what follows, we shall consider two scenarios. The 
first option is to perform the displacement-enhanced pho- 
ton subtraction on one of the two modes only. In this 
case, the filtering operation reads 



Fi = D\{a)aDA{a) ®Ib = {a + a) (gjli 



(3) 



where I denotes the identity operator. Secondly, we can 
subtract photons from both modes, which is experimen- 
tally more challenging, but brings another degree of free- 
dom that can be used to optimize the protocol. In this 
situation, we get the filtering operator 



- D\{a)aDA{a)®D\{p)bDB{P) 
= (a + a) (g) {h + P). 



(4) 



III. WEAK INITIAL SQUEEZING 

In this Section, we assume that the input state is only 
weakly squeezed, A <C 1, so that it can be approximated 
as 
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(t2|2,0)+r2|0,2)). (5) 



This truncated expansion contains the dominant vacuum 
term and a two-photon contribution originating from the 
two-photon component of \s) split on the beam split- 
ter BS. To quantify the entanglement of the filtered 




FIG. 1: Entanglement concentration by photon subtraction 
enhanced by local displacements. The input state is obtained 
by mixing a single-mode squeezed vacuum state with vacuum 
on beam splitter BS. Prior to photon subtraction, local dis- 
placement Da{oi) [DbIP)] is applied, which is undone after 
the subtraction by D\{a) [_D^ (/?)]. The subtraction is ac- 
complished using a highly unbalanced beam splitter with an 
auxiliary mode C (D) in the vacuum state. A click of the 
on-off detector heralds a successful photon subtraction. The 
photon subtraction can be performed either only on one, or 
on both of the modes A, B. 
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FIG. 2: Dependence of entropy of entanglement Es of trans- 
formed state on coherent displacement a with single-mode fil- 
tration. Input states were obtained by mixing the SMSV state 
with vacuum on a beam splitter with reflectance (a) r = l/\/2 
(b) and r = l/\/3. The input squeezing is A = 0.01. 



state \iPab) <x. Fji/'in), we use the entropy of entangle- 
ment Es defined as the von Neumann entropy of one 
of the reduced states, Es = S{pa) = S{pb), where 
PA = T^t:b 114' ab){^ab\], and pB is defined similarly. 
For the single- mode filtration ([3|), the resulting state 
= AlV'in) reads 



1^ 



1/ ~ a|0,0) + Art|0, 1) + At^|l,0) 
aX 



+ aXrt\l,l) + ^(i2|2,0) +r2|0,2)). (6) 



This formula suggests that the best choice of displace- 
ment is a = 0, together with using a balanced beam 
splitter for the initial state preparation, as this gives a 
maximally entangled two-qubit state -^(|0,1) -I- |1,0)), 
i.e., a single-photon state evenly split between modes A 
and B. We have carried out extensive numerical calcula- 
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tions that support this conjecture. As an example, ob- 
tained by numerical calculations truncated at the Fock 
number rimax = 10, Fig. [2] shows the dependence of 
Es on displacement a for two different beam splitter 
reflectances, namely, r = 1/^/2 (left) and r — 
(right). Independent of the beam splitter reflectance, 
we get maximum entanglement for zero displacement. 
Furthermore, in this case the normalized filtered state 
is given by r|0, 1) + t|l,0), whose entanglement is maxi- 
mized for r = t — l/\/2. 

With photon subtraction performed on both modes, 
we get another degree of freedom, the second coherent 
displacement /3, which can be used to maximize the en- 
tanglement of the output state. Using the input-state 
approximation the state \ip2) = ^2|V'in) reads 

IV'2) « {Xrt + af3)\0,0)+X{ar + f3t){t\l,0)+r\0,l)) 

+ ^aPit^\2,0) + V2rt\l,l) +r^\0,2)). (7) 



We can see that the entanglement is maximized for dis- 
placements near the hyperbole a/? = —Art. If this condi- 
tion is satisfied, the vacuum term in state |'02) disappears 
due to destructive interference. For A ^ 1 the state then 
effectively becomes the entangled single-photon state, 
\ip2) ~ i|l,0) + ^10,1) because the last term becomes 
proportional to A^, hence negligible. This may look dis- 
appointing because this state can be obtained much more 
easily by subtracting only a single photon. 



However, there is a further narrow peak clearly visible 
in Fig. |4] which emerges at a = ^/Xt and /3 = — \/Ar. At 
this point, the single-photon part of \1jj2) is also elimi- 
nated by destructive interference as ar -I- /3t = 0, so the 
two-photon part of 1^2) becomes dominant. Since this 
part is proportional to A^ , we must take into account the 
four-photon terms in the expansion of \ipin) in our analy- 
sis. After photon subtractions, the dominant part of this 



The entanglement of the output state is shown in Fig. [S] term which is proportional to A^ reads 
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|4-fc,fc) = -^X\tit^\2,0) +V2rt\l,l) +r^\0,2)). 
v2 



(8) 



r 



Remarkably, this term has the same structure as the two- 
photon part of IV'2) in Eq. ©• Therefore, if a = y/Xt and 
/3 = — -s/Ar the state after filtration becomes an entangled 
two-photon state 

1^2) ~ V2rt\l,l)+t^\2,0)+r^\0,2). (9) 

The dependence of the entropy of entanglement of this 
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FIG. 3: The entropy of entanglement for the two-mode sub- 
traction scheme. The output states with increased entangle- 
ment lie near a hyperbole = —Xrt/a. On the hyperbole, the 
dominant vacuum term vanishes completely, giving maximum 
output entanglement. The input squeezing is A = 0.01 and 
the use of a balanced beam splitter, r = l/-\/2, is assumed. 



state on the beam splitter refiectance is plotted in 
Fig-IHb). For a balanced beam splitter, we get Es = 1.5, 
which is not far from the maximally entangled qutrit 
state, for which we have Es = log2 3 = 1.585. We can 
conclude that the coherent displacements improve the 
performance of the scheme with two subtractions and al- 
low us to generate more entanglement than the scheme 
with a single subtraction only. The price for this entan- 
glement gain is the reduced success rate that scales as A^ 
for single subtraction while it becomes proportional to 
A^ for double subtraction. The success probability and 
performance of the protocol are further analyzed in the 
next section where we extend our investigation beyond 
the weak-squeezing limit. 



IV. ENTANGLEMENT CONCENTRATION 
WITH ARBITRARILY SQUEEZED STATES 

We now provide a more generic analysis allowing for 
arbitrarily strong squeezing and accounting for the ac- 
tual implementation of photon subtraction as sketched 
in Fig. [1] We shall assume that the tap-off beam split- 
ters have amplitude transmittance ts and reflectance rs 
with + r| = 1. The avalanche photodiodes commonly 
employed in this kind of experiments cannot distinguish 
the number of detected photons and respond with a bi- 
nary outcome (click or no click). Consequently, the out- 
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FIG. 4: Two-mode entanglement concentration with recipro- 
cal displacements, /3 = —Xrt/a. In panel (a), the dependence 
of the entropy of entanglement on displacement a is plotted 
for a balanced beam splitter. The plateaus with Es ~ 1 cor- 
respond to the qubit state t|l,0)-|-r|0,l), while a narrow peak 
appears for the qutrit state (O when the one photon terms 
in (O vanish due to destructive interference. The entropy of 
entanglement of the qutrit state (|9]) as a function of the beam 
splitter intensity reflectance 7i = is represented in panel 
(b). 



with ko = inax(0, m — n) is the matrix element of the dis- 
placement operator D{a). If the detector which heralds 
successful photon subtraction only distinguishes between 
presence and absence of photons, the resulting state of 
modes A and B is given by a mixture of all states corre- 
sponding to subtraction of one or more photons. 
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(15) 
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The states are not normalized and the norm of 

I'f/''^'^^) is equal to the probability of subtracting k pho- 
tons. The overall success probability therefore reads 

A limited efficiency rj of the heralding detectors can be 
modeled by a beam splitter with intensity transmittance 
rj followed by a perfect detector with unit efficiency. If k 
photons are reflected on the tap-off beam splitter, then 
the inefficient detector does not register any photon with 
probability (1 — r;)*"', and clicks with probability 1 — (1 — 
1])''. The output state thus reads 



put state after photon subtraction becomes a mixture of 
states corresponding to subtractions of different numbers 
of photons. Since our goal is to evaluate the entangle- 
ment of the output state, we continue the calculations in 
the Fock state representation. We first derive the output 
state for single-mode photon subtraction; the extension 
to the two-mode subtraction is straightforward. 

It is convenient to rewrite the input two-mode state as 
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where 



0, 



m + n even, 



(^)! V^.' 

m + n odd. 
(11) 

Let us now determine a pure state of modes A 

and B after coherent displacement DA{a) and conditional 
subtraction of k photons from mode A. After some alge- 
bra one finds that 



m.n— 



where 
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(12) 



(13) 



(14) 



Piiv) = 



(16) 



k=l 



where Pi (77) = EZii^ " (1 " v)'']{'4'^''^\i^'^''^) ■ 

Going beyond transformation of specific input states, 
we next provide an explicit expression for a non-unitary 
quantum filter describing the combined effect of coherent 
displacement and subtraction of exactly k photons from 
the mode. Similarly as in preceding sections, it is help- 
ful to assume that the coherent displacement D{a) is 
undone after subtraction. Taking into account the re- 
fiectance of the tap-off beam splitter, the optimal in- 
verse displacement reads D{—tsa). The formula for 
the filtering operator P*^'^' is most easily derived by 
considering an input coherent state I7). We use the 
properties of displacement operator I)(7)|0) — I7) and 
D{a)D{'j) = e("'^*-"*'')/2z)(a -1-7) and the transforma- 
tion of the coherent state on the tap-off beam splitter, 
\Oa\0)c \tsOA\rsOc, where C is the auxiliary vac- 
uum input port of that beam splitter. In this way we 
obtain 

F«|7) - e4("T*-"*^)c(fc|rs(a + 7))c|ts7)A 



g-#-("7*-a'7)g- 



Vki 



which is valid for arbitrary I7). From this expression we 
can determine the filtering operator. 



'(a + a)* 



(17) 



If we set a = we recover filter describing fc-photon 
subtraction 33] and in the limit — > we obtain the kih 
power of displaced annihilation operator (a -I- a)''. The 
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factor describes state attenuation on a beam splitter 
with transmittance ts and an additional term e~''s" " 
emerges for nonzero displacement a. 

The output state for the two-mode photon subtraction 
can be determined in a similar manner as for the single- 
mode subtraction. After some algebra, we get 



P2 



■'^ " k,i=i 



(18) 



where state = X^m n=o 1"^' obtained 

after subtracting k photons from mode A and I photons 
from mode B, 
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FIG. 5: Performance of entanglement concentration for the 
single-mode subtraction scheme. Panel (a) shows logarithmic 
negativity as a function of displacement a for initial squeez- 
ing A = 0.2. Zero displacement gives, as in the case of weak 
squeezing, the most entanglement on the output. Panels (b) 
and (c) show the logarithmic negativity and success probabil- 
ity, respectively, as functions of the initial squeezing A for zero 
displacement. The results were obtained for Rs — r% — 0.1, 
considered detection efficiencies are rj — 1 (solid line), rj = 0.5 
(dashed) and rj = 0.1 (dotted). 



a,fc=0 



and 



k,l=l 

Note that we assume the same reflectance rg of both tap- 
off beam splitters and equal detection efficiency rj of both 
detectors. 

To quantify the entanglement of the output mixed 
state, we use logarithmic negativity, defined as [s^ 



En{p) = log2 



Ip^-lll, 



(20) 



which is a computable entanglement measure for mixed 
states. Here, Ta denotes partial transpose of density ma- 
trix p with respect to mode A and || • ||i denotes trace 

norm, ||X||i = TtVx'^X. 

In numerical calculations, the infinite summations are 
truncated at sufficiently large Fock number rimax, where 
J^max = 10 is sufficient for the values of squeezing consid- 
ered in this paper. Numerical results obtained for pho- 
ton subtraction performed on one mode are plotted in 
Fig. m We consider balanced beam splitter for prepa- 
ration of the input entangled state. In fact, the out- 
put entanglement is maximized for a slightly unbalanced 
beam splitter, as the single-mode filtration is asymmet- 
ric, but the change of reflectance is negligible (approx- 
imately 0.5 % higher intensity reflectance for A — 0.2 
and r| = 0.1). As in the case of weakly squeezed in- 
put states, the most entanglement is achieved with zero 
displacement; see Fig. [Sja). In the weak-squeezing limit. 
En is limited by the reflectance of the tap-off beam split- 
ter and the detector efficiency and more than 0.9 e-bit of 
entanglement can be extracted even for realistic param- 
eters r| = 0.1 and r] = 0.1. As shown in Fig. [SJb), the 
maximum entanglement grows only weakly with increas- 
ing squeezing (A = 0.4 corresponds to more than 3 dB 
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FIG. 6: Entanglement concentration with photon subtraction 
on both modes. Contour plot (a) shows dependence of the 
logarithmic negativity on the local displacements a, (3, while 
in panel (b), the cut along /3 = —a is plotted for positive 
values of a. The parameters used are A = 0.2, Rs = 0.1, 

77=1. 



of squeezing). On the other hand, the amount of squeez- 
ing has a strong impact on the success probability of the 
photon subtraction. In the considered setup with highly 
a unbalanced tap-off beam splitter, the limited efficiency 
of the trigger detector mainly decreases the success prob- 
ability that is approximately proportional to rj, while the 
entanglement is not much reduced. 

The performance of entanglement concentration with 
photon subtraction on both modes is plotted in Fig. [HI 
The resulting dependence of entanglement on the coher- 
ent displacements a and /3 is similar to that of weakly 
squeezed states (cf. Fig. [3]). In what follows, we restrict 
ourselves to the diagonal line /3 = —a [Fig. [Hb)], at 
which the maximum value of logarithmic negativity lies. 
Numerical simulations also indicate that the restriction 
to real displacements is not limiting, as complex phase 
factors do not increase the entanglement. In Fig. [71 we 
investigate optimal entanglement concentration for dif- 
ferent values of initial squeezing A. The figure shows 
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FIG. 7: Optimal two-mode filtration in dependence on the 
squeezing constant A. Panel (a) shows the optimal displace- 
ment Oopt, panel (b) represents the logarithmic negativity 
En, and in panel (c), we plot the success probability P2- The 
results were obtained for Rs = 0.1, and rj — 1 (solid line), 
77 = 0.5 (dashed), rj = 0.1 (dotted). The optimal displace- 
ment is mostly unaffected by limited detector efficiency; for 
better readability, only line for = 1 is shown in panel (a). 
The dot-dashed line in this panel represents the \/X/2 ref- 
erence curve, the optimal displacement for weakly squeezed 
input states. 

optimal displacement aopt, the resulting logarithmic neg- 
ativity En and the success probability P2. Similarly to 
the single-mode subtraction, the initial squeezing influ- 
ences mainly the success probability while the logarith- 
mic negativity only weakly increases with A. Compared 
to the single-mode subtraction scheme, the amount of 
entanglement we can get from the system is higher, at 
a cost of smaller success probability. Inefficient detec- 
tion lessens the entanglement of the output state more 
than in the case of the single-mode subtraction, as two 
detectors are used. As shown in Fig. El^a), for A < 0.2 
the optimal displacement aopt practically coincides with 
the optimal value a/A/2 determined analytically in the 
weak-squeezing limit, while for larger A the optimal dis- 
placement exceeds this analytical prediction. The success 
probability is the parameter that is affected most by im- 
perfect detection and approximately scales as P2 oc rj^. 
Obviously, this is the main disadvantage of the two-mode 
subtraction compared to the single-mode subtraction. 

V. CONCLUSIONS 

In summary, we have analyzed continuous variable en- 
tanglement concentration based on photon subtraction 



enhanced by local displacements. We considered split 
single-mode squeezed vacuum state as input of the pro- 
tocol, which makes the whole setup easier to implement 
experimentally and directly corresponds to the experi- 
ment by Takahashi et al. [1^. We have investigated 
and compared two different scenarios - in the first one, 
the displacements and photon subtraction are performed 
on one mode only, which makes it experimentally easier, 
while the other scheme employs the filtration operation 
on both modes and thus introduces another degree of 
freedom that can be used to extract more entanglement 
in case of successful photon subtraction. 

In the case of the single-mode photon subtraction lo- 
cal displacements do not provide any improvement. From 
this point of view, simple photon subtraction, as already 
demonstrated by Takahashi et al. [2^, is the best op- 
tion, although other types of local operations, such as 
squeezing, might work differently. On the other hand, the 
scheme with photon subtraction on both modes can be 
significantly improved by local displacements. Moreover, 
the entanglement of the resulting state can be higher than 
in the case of the single-mode subtraction. This is a re- 
markable improvement, since single-photon subtraction 
can outperform photon subtraction on both modes with- 
out displacements. 

The required combination of photon subtraction and 
coherent displacement has already been successfully 
demonstrated experimentally and explored for engineer- 
ing of arbitrary qubit states formed by superpositions of a 
squeezed vacuum and squeezed single-photon states [s^ . 
In this context, we note that the entanglement concen- 
tration scheme studied in this paper not only enhances 
the entanglement of the state but also modifies its struc- 
ture. Thus, the considered protocol can also be seen as 
an example of a bipartite quantum state engineering by 
local operations. In the future, it would be interesting 
to investigate the power and limitations of this approach 
and determine which states can be prepared from initial 
Gaussian entangled states by local photon subtractions, 
local Gaussian operations and classical communication. 
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